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Abstract. From 43 GHz VLB A observations of the pair 
of radio sources QSO 1928+738 and BL 2007+777 we have 
demonstrated the feasibility of precision phase-delay dif- 
ferential astrometric techniques at millimeter wavelengths. 
For a pair of sources with 5° separation and high antenna 
elevations, we have shown that present astrometric mod- 
els and millimeter arrays are advanced enough to model 
the differential phase-delay to within 2 picoseconds, less 
than one tenth of a phase-cycle at 43 GHz. The root-mean- 
square of the differential phase-delay residuals is domi- 
nated by the fluctuations of the atmospheric water vapor. 
We have determined the relative position of the observed 
sources with a precision twofold better than previous de- 
terminations at lower frequencies and, more importantly, 
largely free from ambiguous definitions of the reference 
point on the structure of the radio sources. Our result 
makes 43 GHz VLBI phase-delay differential astrometry 
an ideal tool to study the absolute kinematics of the highly 
variable structures of regions near the core of extragalactic 
radio sources. 

Key words: astrometry - techniques: interferometric - 
quasars: individual: QSO 1928+738 - BL Lacertae object: 
individual: BL 2007+777 



1. Introduction 

One of the trends in Very-Long-Baseline Interferometry 
(VLBI) is to augment the angular resolution of the ob- 
servations in search for a more detailed view of the inner 
structure of extragalactic radio sources. This is effectively 
carried out by either observing at millimeter wavelengths 
(mm- VLBI) or, at cm- wavelengths, by combining ground 
telescopes with antennas in space. The correct interpre- 
tation of these high-resolution observations is of much 
relevance since they map the morphology of highly vari- 
able regions close to the central engine of AGNs. How- 
ever, multi-epoch analyses directed to understand the dy- 
namical behavior of these inner regions critically depend 



Send offprint requests to: J.C. Guirado 



on the alignment of the images: no solid conclusions can 
be extracted without an accurate source component (i.e. 
core) identification. In particular, VLBI reveals that cm- 
wavelength components break up in complex structures 
with multiple features at mm- wavelengths. These com- 
pact mm-features show a strong variability, which may 
be the result of phenomena only seen so far in numerical 
simulations (Gomez et al. 1995). For a meaningful phys- 
ical understanding of those compact features, a detailed 
knowledge of the (absolute) kinematics of the region is 
crucial. It is therefore highly desirable to extend precision 
differential phase-delay astrometry to mm-wavelengths. 
In this Letter we demonstrate the feasibility of us- 
ing phase-delay differential astrometry at 43 GHz. We 
have selected the pair of sources QSO 1928+738 and 
BL 2007+777, —5° apart, with flat spectra, high flux den- 
sities, and rich structures at 43 GHz. The astrometry anal- 
ysis of these data show the advantages and possibilities of 
mm-wavelength differential astrometry. 



2. Observations and Maps 

We observed the radio sources QSO 1928+738 and 
BL 2007+777 at 43 GHz on 1999 January 3 from 15:00 
to 23:30 UT. We used the complete Very Long Base- 
line Array (VLBA) recording in mode 256-8-2 in left cir- 
cular polarization to achieve a recording bandwidth of 
64 MHz. We interleaved observations of QSO 1928+738 
and BL 2007+777 using integration times between 40 and 
130 seconds on each source to make total cycle time dura- 
tions between 110 and 300 seconds (antenna slew time was 
15 seconds). The data were correlated at the National Ra- 
dio Astronomy Observatory (NRAO, Socorro, NM, USA). 
Detections were found on both sources to all stations but 
Hancock, presumably due to severe weather conditions 
at this site. We made manual phase calibration, visibil- 
ity amplitude calibration (using system temperatures and 
gain curves from each antenna), and fringe fitting at the 
correlation position with the NRAO Astronomical Image 
Processing System (AIPS). For astrometric purposes, we 
further processed the data in AIPS (tasks MBDLY and 
Cl2hf) to obtain, for each baseline and epoch, estimates 
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Fig. 1. 43 GHz hybrid maps of QSO 1928+738 (left) and 
BL 2007+777 (right) at epoch 1999.01. Contours are - 
0.5,0.5,1,2,4,8,16,32,64, and 90% of the peak of brightness, 
0.56 Jy/beam, for QSO 1928+738, and -1,1,2,4,8,16,32,64, 
and 90% of the peak of brightness, 0.22 Jy/beam, for 
BL 2007+777. The restoring beam (shown at the bot- 
tom left corner of each map) is an elliptical Gaus- 
sian of 0.24x0. 14mas (P.A. 21°) for QSO 1928+738, and 
0.23x0.14mas (P.A. 29°) for BL 2007+777. The astromet- 
ric reference points are labeled as Ql, and correspond to 
the peak of brightness of each map. See text. 

of the group delay, phase-delay rate, and fringe phase at 
a reference frequency of 43,185 MHz. We discarded data 
from Saint Croix and Pie Town stations, which showed 
unacceptable scatter in the observables, for the astromet- 
ric analysis presented in Sect. 3. For mapping purposes, 
we transferred the data into the Caltech imaging program 
difmap (Shepherd et al. 1995). We performed several it- 
erations of self-calibration in phase and gain. We present 
the resulting hybrid maps in Fig. 1. 

The 43 GHz map of QSO 1928+738 displays several jet 
components extending southwards. All these features ap- 
pear blended together as only one or two components in 
previous maps at cm-wavelengths (Guirado et al. 1995, 
hereafter G95; Ros et al. 1999, hereafter R99), and those 
obtained with space VLBI (Murphy et al. 1999). The 
brightest knot, labeled as Ql, is probably a jet compo- 
nent, unless the source is two-sided, but it is compact and 



well defined. Thus, it constitutes an appropriate reference 
point for relative astrometry at a single epoch. However, 
this component would not be a suitable reference point 
for a multi-epoch comparison of the relative separation 
between the two sources as it is likely to move and evolve 
in brightness and shape over time. 

The 43 GHz map of BL 2007+777 represents a significant 
improvement in the knowledge of the inner structure of 
this source (see Fig. 1). The brightest knot seen earlier at 
cm- wavelengths (Guirado et al. 1998, hereafter G98; R99) 
breaks-up in at least three new features. The kinematic na- 
ture of component Q2, almost as bright as the easternmost 
component Ql, is of much interest. The brightest feature 
of BL 2007+777 at cm- wavelengths, a blending of all com- 
ponents seen within 1 mas from the origin in our 43 GHz 
map, has been taken as a reference point for astrometry 
(G95; G98; R99); even more, this feature has been consid- 
ered stationary for multi-epoch astrometry analyses. Ac- 
cordingly, should component Q2 be a travelling knot, the 
selected reference point at cm-wavelengths is likely to be 
not stable over time and part of previous astrometry re- 
sults must be revised. 



3. Astrometry Analysis 

A goal of this research has been to calibrate the limitations 
of our standard astrometry procedure for 43 GHz, as well 
as to study the potential precision of the astrometric data 
at this frequency. Therefore, the data-reduction procedure 
for the 43 GHz observation deliberately followed the same 
steps as those used for the 5 GHz (G95) or 8.4 GHz obser- 
vations (G98; R99). We briefly go again over each step of 
this analysis: For our 43 GHz data, (i) we predicted, via a 
precise theoretical model of the geometry of the array and 
the propagation medium, the number of cycles of phase 
between consecutive observations of the same source to 
permit us to "connect" the phase delay (e.g. Shapiro et 
al. 1979; G95; R99); (ii) we defined as reference points in 
the 43 GHz images of the two radio sources the maximum 
of the brightness distribution (components Ql in the maps 
of QSO 1928+738 and BL 2007+777; see Fig. 1) and sub- 
tract the contribution of the structure of the radiosources, 
with respect to the reference points selected, from the 
phase delays; (iii) we formed the differenced phase de- 
lays by subtracting the residual (observed minus theoreti- 
cal values) phase delay of BL 2007+777 from the previous 
observation of QSO 1928+738; and (iv) we estimated the 
relative position of the reference points from a weighted- 
least-squares analysis of the differenced residual phase de- 
lays. For this analysis we used an improved version of the 
program VLBI3 (Robertson 1975). 

In step (i), the geometry of our theoretical model (set of 
antenna coordinates, Earth-orientation parameters, and 
source coordinates) was consistently taken from IERS 
(IERS 1998 Annual Report, 1999). The theoretical model 
also accounted for the effect of the propagation medium 
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Table 1. Contributions to the standard errors of the es- 
timates of the coordinates of QSO 1928+738 minus those 
of BL 2007+777 ( SAa, 5 AS) from the sensitivity study. 
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a The standard deviation of the fixed geometrical parameters of our 
astromctric model (all entries but ionosphere) were taken from IERS 
Annual Report 1998 (1999). The 2 cm standard deviation of the site 
coordinates corresponds to each of the three coordinates for each 
antenna site. 

b Notice that the values of <5Aa are in /is. To convert fis to fias, the 
factor 15- cos (5qso 1928+738 (~4.2 ) must be used. 

c The statistical errors include the uncertainties of the troposphcric 
zenith delays at the nodes of the piecewise linear function used in the 
troposphere model (see Sect. 2). 

d Standard deviation provided by the global ionospheric maps at each 
site. 1 TECU = lxlO 16 elm" 2 . 

e Root-sum-square of the tabulated values. 



in the astrometric observables. We modeled the iono- 
spheric delay by using total electron contcct (TEC) data 
from GPS-based global ionospheric maps generated at the 
epoch of our observations by the Center for Orbit Deter- 
mination in Europe (Schaer et al. 1998). We followed the 
geometric corrections described in Klobuchar (1975) and 
Ros et al. (2000). We modeled the tropospheric zenith 
delay at each station as a piecewise-linear function char- 
acterized by values specified at epochs one hour apart. We 
calculated a priori values at these nodes from local surface 
temperature, pressure, and humidity, based on the model 
of Saastamoinen (1973). The antenna elevations were al- 
ways higher than 20° at all stations; this allowed us to use 
the dry and wet Chao mapping functions (Chao 1974) to 
determine the troposphcric delay at non-zenith elevations 
for each observation at each site. We estimated the tro- 
posphcric zenith delay at the nodes of each station, along 
with the relative position of the sources, from a weighted- 
least-squares analysis. 

4. Results and Discussion 

From the astrometric analysis described in Sect. 3, we ob- 
tain the following J2000.0 coordinates of QSO 1928+738 
minus those of BL 2007+777 at 43 GHz: 



Aa = -0 h 37 m 42? 503443 ± 0?000026 
AS = -3° 54' 41 "677208 ± 0"000056 

where the quoted uncertainties are overall standard errors 
(see Table 1), nearly twofold smaller than the standard 
errors corresponding to previous determinations at 5 and 
8.4 GHz. From the comparison of the results of the sensi- 
tivity analysis displayed in Table 1 with similar sensitivity 
analyses at lower frequencies (G95; R99), we see that the 
improvement in precision comes from (i) the small contri- 
bution to the standard errors of the reference point iden- 
tification in the map (dominated by image noise), as a 
consequence of the improvement of resolution of the maps 
and of the lack of ambiguity in selecting the components 
acting as reference, and (ii) the negligible contribution of 
the ionosphere, that scales down by a factor of 25 with 
respect to its contribution at 8.4 GHz. As occurs at cm- 
wavelengths, the quoted standard errors of the relative 
position are dominated by the uncertainties of the fixed 
parameters of the astrometric model (entries 3 to 10 of 
Table 1), and, in particular, by the uncertainties of the co- 
ordinates of the reference source, as expected for objects 
with a large angular separation (notice that this error is 
not frequency dependent). The comparison and interpre- 
tation of the relative position estimate at 43 GHz with 
previously reported estimates at lower frequencies will be 
postponed to a later publication where the comparison 
will be made in great detail. 

The postfit residuals of the differenced phase delays cor- 
responding to all baselines included in our analysis are 
shown in Fig. 2. Note the scale of the plots, ±23 ps, cor- 
responding to ±1 phase cycle. The average root-mean- 
square (rms) of the postfit residuals is 2ps, less than one 
tenth of the phase cycle at 43 GHz. At this level of pre- 
cision, the absence of systematic effects validates both 
the astrometric model, based on IERS standards, and 
the propagation medium procedures for mm-wavelength 
VLBI astrometry (at least for cycle times, source separa- 
tions, weather conditions, and antenna elevations similar 
to those presented in this paper) . To calibrate the quality 
of our procedure, we compared the residual of the dif- 
ferenced phase delay with similar residuals corresponding 
to observations at 8.4 GHz and 5 GHz made in the past 
(G95, G98). The rms of the residuals are about 15, 9, and 
2ps for the data sets at 5, 8.4, and 43 GHz, respectively, 
which expressed in equivalent-phase yield postfit residu- 
als of roughly 30 degrees at each of the three frequencies. 
This similarity of the rms expressed in phase at all the ob- 
served frequencies demonstrates not only that the phase 
connection process is feasible at 43 GHz, but which is also 
of no less quality than at lower frequencies. 
Likely, the most important contributors to the scatter of 
the phases at 43 GHz are the unmodeled variations of the 
refractivity of the neutral atmosphere. From the average 
rms of 2 ps of the phase residuals of Fig. 2, and assum- 
ing uncorrelated contributions from the antennas forming 
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each interferometric pair, the average uncertainty for the 
single-site phase delay is \[2 ps. This uncertainty is in good 
agreement with the predictions of water vapor fluctuations 
on time scales of 100 seconds (~2ps) based on refractivity 
patterns described by Kolmogorov turbulences (Treuhaft 
& Lanyi 1987). 

The importance of our result translates to VLBI phase- 
referencing mapping. This technique (see e.g. Lestrade et 
al. 1990) relies completely on the behavior of the phase 
of the reference source (usually a strong radio emitter) to 
predict the phase of the target source (usually a weak radio 
emitter). Beasley & Conway (1995) provide useful expres- 
sions for the maximum cycle time for phase-referencing 
with the VLBA. Under good weather conditions, average 
antenna elevation of 40°, and with the requirement that 
the rms phase between scans is <90°, the maximum cycle 
time at 43 GHz is ^100s. This estimate should be short- 
ened if atmospheric spatial variations from different lines 
of sight are considered. Actually, the facts are more favor- 
able. For sources separated 5° on the sky, high antenna 
elevations, and good weather conditions, our results show 
that (i) the rms of the phases is below 90° throughout 
the experiment and does not seem to be substantially de- 
pendent of the cycle times used during our observation 
(100-300s); and (ii) the expected average uncertainty in 
interpolating the phases of one source to the epoch of the 
other is ~30° in the differenced phase. This value is not 
larger than usual phase errors in phase-reference mapping 
at cm- wavelengths (e.g. Lestrade et al. 1990). Therefore, 
with the proper cycle times and nearby calibrator sources, 
diffraction limited VLBI phase-reference images at 43 GHz 
should be possible. 

Our observations have shown that VLBI differential as- 
trometry at 43 GHz provides high-precision relative posi- 
tions. At this frequency, the astrometric precision is nearly 
equivalent to the resolution of the maps, and the reference 
point selected in the source structure might be associated 
with the core. This makes 43 GHz differential astrometry 
an ideal technique to trace unambiguously the kinematics 
of the inner regions of the extragalactic radiosources. 
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